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We investigate the transmission properties of quasiperiodic or aperiodic structures 
based on graphene arranged according to the Cantor sequence. In particular, we have 
found self-similar behaviour in the transmission spectra, and most importantly, we 
have calculated the scalability of the spectra. To do this, we implement and propose 
scaling rules for each one of the fundamental parameters: generation number, height 
of the barriers and length of the system. With this in mind we have been able to 
reproduce the reference transmission spectrum, applying the appropriate scaling rule, 
by means of the scaled transmission spectrum. These scaling rules are valid for both 
normal and oblique incidence, and as far as we can see the basic ingredients to obtain 
self-similar characteristics are: relativistic Dirac electrons, a self-similar structure and 
the non-conservation of the pseudo-spin. This constitutes a reduction of the number 
of conditions needed to observe self-similarity in graphene-based structures, see Dfaz- 
Guerrero et al. [D. S. Dfaz-Guerrero, L. M. Gaggero-Sager, I. Rodriguez-Vargas, and 
G. G. Naumis, arXiv:1503.03412vl, 2015]. 
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I. INTRODUCTION 


In the last decade, the aperiodic modulation has attracted great interest due to its pe¬ 
culiar electronic and optoelectronic properties^ such as the transmission properties that 
appear highly fragmented displaying self-similar patterns along with a fractal-like electronic 
spectra,- different to what occurs in periodic superlattices. All these characteristics are a 
direct consequence of the quasiperiodicity that underlies in these aperiodic structures. 

Quasiperiodic or aperiodic systems are structures that can be considered as fractals,- ho¬ 
mogeneous and self-similar geometrical objects, that also are defined as structures ordered 
but without periodicity. In this context, the aperiodic structures or also commonly known 
as deterministic media are classified in two group according to the constructions algorithm. 
Those generated by means of a substitutional procedure named as substitutional sequences 
(Fibonacci, Thue-Morse, Double-Period, Rudin-Shapiro and Silver-Mean) and those gener¬ 
ated in a similar way to fractal sets,- named as model fractal structures (Cantor sets and 
Koch fractals). Among the peculiar characteristics of the quasiperiodic structures excels the 
so-called self-similarity.- This property means that one or more parts of a structure or object 
have exactly or approximate the same shape of the whole, and even more if these parts are 
magnified at different scales the result is the entire object. 

Self-similarity is found everywhere, from geographic structures, statistical natural pro¬ 
cesses to complex networks.- In terms of length scales, we can say that self-similarity is found 
in systems of the size of the atom to systems of the size of the universe.- Specifically, this out¬ 
standing property was recognized in the coasts, in the distribution of species, in the growth 
of trees, rivers and lungs, to mention a few.- - - Within this context, the aforementioned self¬ 
similar parts of a structure can be reproduced at different scales, so it is natural to assume 
that there are scaling rules that relates the scalability of thereof. However, in most of the 
multiple reports that we can find in the literature about solid-state systems, self-similarity 
is let as a matter of visual perception, and rarely the scaling rules are reported. 7 

In this sense, since the discovery of graphene^ - — researchers have begun to deal with ape¬ 
riodic structures based on this prominent material.— - — For instance, the Fibonacci sequence 
is studied in monolayer— and bilayer— graphene superlattices. In particular, a striking self¬ 
similar behaviour is reported for different generations number of the aperiodic sequence,— it 
is also found that the transmission spectra are fragmented, and even more appear in groups 
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and the low-energy conductance is considerably amending.— With respect to the Cantor set, 
a few works have been devoted in the same direction.— - — Specifically, a magnetic modula¬ 
tion has been implemented to study the transport properties through graphene-based fractal 
barriers.— Despite the particularities of these reports, we can see that they have a common 
aspect, a quasiperiodic structure with self-similar properties, which are manifested in the 
transmission spectra, ffowever, the scaling rules of the self-similar patterns are not found, 
and once again self-similarity is considered as a matter of visual perception. Thus, a simple 
question arises, what are the appropriate scaling rules that can reproduce the self-similar 
patterns?. 

Following the lines of a previous work,— we obtain the scaling rules of the self-similar 
transmission patterns of Cantor graphene structures (CGSs). In this case CGSs are con¬ 
structed by simply scaling (spatially) the barriers of the system according to the triadic 
Cantor set rules. Scaling rules between generations, lengths of the system and heights of 
the barriers are obtained. These scaling properties are valid for both normal and oblique 
incidence. We also carry out a comparative analysis of the scaling expressions in substrate- 
based CGSs (SCGSs), electrostatic CGSs (ECGSs) and Al x Gai_ 3 ,As/GaAs semiconductor 
structures (SSs), finding a better scaling in the case of SCGSs. Our results indicate that 
conditions like the non-symmetry (spatially speaking) and the scaling in the energy axis of 
the structure are not necessary to obtain self-similarity in the transmittance.— Therefore, 
as far as we can see, the conditions to obtain self-similarity in the transmission properties 
are: relativistic Dirac electrons, a self-similar structure and the non-conservation of the 
pseudo-spin (electrons in SCGSs). 


II. METHODOLOGY 

Prior to proceeding with the mathematical apparatus that we employed, we will present 
a simple procedure, an iterative method, to obtain the triadic Cantor set sequence. For the 
first generation, N = 1, we choose a segment of any length. For the second generation, 
N = 2, we divide it into three equal parts and then eliminate the central one. Finally, if we 
repeat this procedure at the N — th generation the triadic Cantor set is obtained. Therefore, 
with the aforesaid in the previous section it is possible to develop a fractal arrangement of 
barriers according to the substitution rules of this aperiodic sequence, see Fig. |U [2] and [3l 
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As our main objective is fractal structures in graphene, we will present the fundamentals 
of those structures in which the potential barriers are generated by substrates (Fig. [T]) 
and electrostatic fields (Fig. [2]). The fundamentals of Schrodinger electrons can be found 
elsewhere,— and some minor details can be seen Fig. [3J 

To this regard, in graphene substrates such as silicon carbide (SiC) play an important 
role, since they are capable to change the form of the dispersion relation, from linear to 
parabolic,— as well as to open a bandgap in the energy spectrum.— The Hamiltonian that 
describes this kind of systems is given by, 


H = v F (cr • p) + t'a z , (1) 

where v F is the Fermi velocity of quasi-particles in graphene of the order of c/300, t' = mv 2 F 
is the mass term, cr = (a x , a y ) are the Pauli matrices, p = ( p x ,P y ) is the momentum operator 
and <7 Z the ^-component of the Pauli matrix. This equation can be solved straightforwardly 
giving the following dispersion relation 


E = \Jh 2 Vpq 2 + P 2 , (2) 

here q is the wave vector associated with the SiC region, and t' is proportional to the 
bandgap, E g = 2 t'. Henceforth P=0.13 eV, because we deal with a monolayer graphene 
system. In the case of the wavefunction we have,— 


with 





e ±iqxX+iq v y 


V± 


hv F (±q x + iq y ) 
E + t' 


(3) 

(4) 


where the two components of the wavefunction are referred for the two equivalent triangular 
sublattices A and B in the graphene hexagonal structure, for this reason the wavefunction 
is considered as a bispinor. 

On the other hand, systems with electrostatic barriers can be described by the following 
Hamiltonian, 


H = v F (cr • p) + V(x), 


(5) 
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with V(x) = Vo the one-dimensional potential along the x-direction. Considering the above 
Hamiltonian, the Dirac equation can be solved directly giving the linear dispersion relation, 

E = ±hv F q + Vo, (6) 

where “±” represent electrons and holes, respectively. The wavefunctions are represented 
with the same mathematical form that in the case of substrate barriers, eq. (J3]), but now 
the bispinor coefficients are written as, 


v± = 


hv F (±q x + iq y ) 


(7) 


E-V 0 

Once the explicit form of the dispersion relations and wavefunctions are obtained for both 
substrate and electrostatic structures, it is easy to compute the transmission properties of 
our aperiodic fractal system by means of the transfer matrix formalism. 33,34 Specifically, the 
transmission probability or transmittance conies as, 


T = 


An+ i 

2 1 

Ao 

\M U \ 


( 8 ) 


where the transmitted wave amplitude A ^ + 1 can be calculated in terms of the incident wave 
amplitude A 0 via the transfer matrix, 



= M 

lw+1 i 

UJ 


0 ) 


where the transfer matrix is given by, 


35-37 


(9) 


m = d y 


N 

II OiPjD] 

3 = 1 


-1 


Dn 


( 10 ) 


with Dj and Pj the dynamical and propagations matrices of the structures.— — 


III. RESULTS AND DISCUSSION 

First of all, it is important to emphasize that the aim of our work is to find out the 
scaling properties in the transmission spectra for three fundamental parameters: generation 
N, height of the barriers Vo and length of the system w. So, we study the self-similarity and 
scalability of the transmission spectra of SCGSs (Fig. [T]) , ECGSs (Fig. [2]) and a typical SSs 
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(Fig. [3]). Hereafter, we will focus our attention in three special mathematical expressions 
that can describe the scalability of our systems. 

Before continuing with the presentation of our numerical results, we present important 
aspects of the systems of interest: SCGSs, ECGSs and SSs. In the case of SCGSs, the 
charge carriers are treated as relativistic particles known as massive Dirac electrons. Other 
characteristic is that the Klein effect is ruled out due to the non-conservation of the pseudo¬ 
spin. In the case of ECGSs, the charge carriers are massless Dirac electrons. In this case, 
at normal incidence the Klein effect, perfect transmission, is manifested. At this point, it is 
important to mention that our results for massless Dirac electrons at normal incidence are 
not illustrative, in view of the fact that the transmittance is always T — 1, irrespective of 
the generation, height and length of the system, and the scaling rules will be always valid. 
So, the results for massless Dirac electrons at normal incidence will be not presented. In the 
case of SSs, the charge carriers are Schrodinger electrons. An important characteristic in this 
case is that the transversal motion of Schrodinger electrons is separable, so the transmission 
properties are not dependent of the angle of incidence. Therefore, we will only present results 
for Schrodinger electrons at normal incidence. 

In order to analyse the scaling rules of the transmission spectra between generations, we 
propose an expression that depends on the factor three, which is involved when one changes 
the generation in the triadic Cantor set. Specifically, the expression comes as, 

T n (E) « [T n+1 (E)] 3 . (11) 

Now, we implement this expression for massive Dirac electrons and Schrodinger electrons 
at normal incidence, see Fig. [4] The structural parameters for both systems are the same, 
that is, the height and width of the barriers as well as the generation are the same. In 
particular, we consider barriers of height and starting width of Vo = 0.13 eV and w = 1000a, 
with a = 1.42 A the carbon-carbon distance in graphene, and generations N = 6 and N = 7. 
From now on, N = 6 will be the generation of reference and N = 7 the generation to be scaled 
or simply the scaled generation. In Fig. H](a) we show the transmittance as a function of the 
energy for massive Dirac electrons. The solid-black and the dashed-red curves correspond 
to generations 6 and 7, respectively. As we can see the envelopes of these curves are quite 
similar, the only difference that we can notice is a shift in the vertical axis. Appealing to 
this similarity we implement our scale rule (eq. fTT]) to obtain the transmission spectrum of 
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generation N = 6 by means of the scaling of generation N — 7, that is, Tq(E) m [T 7 (E)] 3 . 
As we can see in Fig. HJF) the matching between the generation of reference (.N = 6) and 
the scaled one (N = 7) is quite good, compare the solid-black and dotted-blue curves. In 
the case of Schrodinger electrons the transmission spectra between generations 6 and 7 are 
also similar, see the solid-black and dashed-red curves in Fig. Q|c). However, in this case, 
in addition to the shift in the vertical axis, we can see an offset between the curves in the 
energy axis. And it is precisely this difference that makes that our scale rule not work well 
in the case of Schrodinger electrons, compare the solid-black and dotted-blue curves in Fig. 
Hd). 

The scaling properties of massive and massless Dirac electrons at oblique incidence (6 = 
7t/4) are presented in Fig. 0 These cases are in contraposition, since in one case the 
conservation of the pseudo-spin is ruled out (massive Dirac electrons), while in the other 
case it is sustained (massless Dirac electrons). The structural parameters are the same as 
in Fig. [4j For massive Dirac electrons, once again, we can see that the transmission spectra 
between generations 6 and 7 are quite similar, the main difference is a relative shift in the 
transmission axis, see Fig. |5](a). By applying our scale rule to T 7 (E), we obtain a very good 
correspondence with respect to the transmission spectrum of reference (Tq(E)), compare the 
solid-black and dotted-blue curves in Fig. [5](b). The main differences between these curves 
are localized in the minimums of the transmittance. In the case of massless Dirac electrons, 
we can notice that the scaling properties are worsen in comparison with the massive case, 
see Fig. [5](d). Even, this case is worse with respect to the the case of Schrodinger electrons, 
compare Figs. QJd) and[5](d). These significant differences between the reference spectrum 
and the scaled one come from the large shifts in both the transmittance and energy axes 
between T 6 (E) and T 7 (E), see Fig. [5](c). 

Continuing with our study about CGSs, we have also found that self-similar patterns 
are visible when two systems that have a different height of the barriers are considered. 
Therefore, the scalability for SCGSs, ECGSs and SSs is examined. Thus, comes the turn to 
view how the heights of the barriers play a fundamental role with respect to the self-similar 
properties of our aperiodic structures. In this regard, the following expression is proposed, 
in general, to analyse the scaling properties between systems with different heights of the 
barriers or simply the scaling between heights,— 
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T Vo (E) « [T. l4 (B)] t, I (12) 

where the subindex Vo indicates the height of the barriers, and k is a factor that comes in 
multiples of two. 

To this respect, Fig. 12 depicts the scaling between heights at normal incidence (6 = 0) 
for both massive Dirac electrons and Schrodinger electrons. In this case the structural 
parameters are N — 6 and w = 1000a. At first, we show comparatively the transmis¬ 
sion spectra for a structure with two different heights of the barriers. In particular, Figs. 
Ha) and (c) correspond to massive Dirac electrons and Schrodinger electrons, respectively. 
Specifically, the transmission spectrum is computed for a height equal to Vo = 0.10 eV and 
Vo = 0.05 eV. Both spectra denoted by T 0 .io(F') and T om (E), solid-black and dashed-red 
curves, respectively. Then, if we implement the scale rule (eq. fT21) for these heights, that 
is, To.io(-E') ~ [T 0 .o 5 (£)] 4 , a scaled transmission spectrum for SCGSs (Fig. Hb)) and for SSs 
(Fig. Hd)) is obtained. If we paid attention of the results for massive Dirac electrons we can 
see a great matching between the reference transmission and the scaled transmission. Thus, 
we can corroborate that the self-similarity and scalability are exhibited in SCGSs, and even 
more interesting, we have found a scaling rule that describe these properties. On the other 
hand, Fig. HH) displays the scaling between heights for the case of Schrodinger electrons. 
The height of the barriers for this system, Al x Gai_ x As/GaAs, is modulated changing the 
aluminum concentration (x). Here, we can see that the reference transmission and the scaled 
transmission, solid-black and dotted-blue curves, have dissimilar envelopes, that is, the scale 
rule is not fulfill at all. Even more, this dissimilarity is easy to visualize in the mentioned 
figure because the reference transmission and the scaled transmission are shifted one respect 
to other in both axes. 

In Fig. [7] we show the scaling properties between heights of the barriers, at oblique 
incidence (9 = 7r/4), for massive and massless Dirac electrons. All the parameters for 
SCGSs and ECGSs are the same as in Fig. [6] The corresponding transmission spectrum for 
massive and massless Dirac electrons is calculated for two different heights of the barrier, 
Figs. Ha) and (c), respectively. At first sight, the transmission spectra for massive Dirac 
electrons seem different, however, they have something in common, that is, both curves 
show a similar envelope, compare the solid-black and dashed-red curves in Fig. Ha). So, it 
is natural to think in a scale rule that reproduces one curve with respect to the other. By 
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applying eq. [12] to the transmittance of height Vo = 0.05 eV, we obtain the corresponding 
scaled curve, dotted-blue curve in Fig. [T^b) . As we can corroborate, a very good matching 
between the reference and scaled transmittance is obtained. The main differences are found 
in the minimums of the transmission. Even, it seems that the matching is better than in 
the case of the scaling between generations, see Fig. E^b). In contrast, Fig. [T^d) clearly 
shows that this scale rule is not suitable for massless Dirac electrons. We can notice a 
remarkable offset between the curves, in both the transmission and energy axis. Even more, 
the envelopes of the curves are not corresponding, see Figs. EK c) and (d). In short, we can 
say that the scaling between heights is valid for massive Dirac electrons, irrespective if the 
electrons impinge at normal or oblique incidence. 

Another self-similar patterns in the transmission properties occur when a structure with 
two different lengths is taken into account. With the length of the system we mean to 
the width of the barrier in the first generation of the structure. Within this context, the 
following expression is used to study this kind scaling property,— 

r„(£) « \T lw (-E)f, (13) 

“ a 

where the subindex w refers to the length of the system, and the factor a comes in multiples 
of two. At this point, it is important to remark that this expression is a little bit different 
with respect to the former ones, since in addition to scale the transmission, it is required 
that the argument (the energy) of it, be scaled as well. 

By applying this expression to the case of massive Dirac electrons and Schrodinger elec¬ 
trons at normal incidence (Fig. |HJ), we can see that a very good scaling is obtained for 
massive Dirac electrons, see Fig. |8](b). In this case the structural parameters are N = 6 and 
Vo = 0.13 eV. We have chosen as lengths of the system w = 1000a and w = 500a, the former 
representing the reference system and the latter the scaled one. In addition, the factor a is 
equal to 2, so the reference and scaled transmissions are related by T W oo a (E) ~ [T 50 oa{^ E)] 4 . 
The matching between these curves is quite surprising, since at first sight they are totally 
dissimilar, compare solid-black and dashed-red curves in Fig. [81(a) . Now considering the case 
for Schrodinger electrons, see Fig. |8](d), turns out that the reference and scaled transmission 
are practically different. Even, this case is worsen as compare with the other scale rules, 
see Figs. |4](d) and [61(d). So, as far as we can see the transmission properties of Schrodinger 
electrons do not manifest a self-similar behaviour. On the other hand, if we apply this 
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scale rule to the propagation of electrons in SCGSs and ECGSs, we find that a very good 
scaling is presented in the case of massive Dirac electrons, see Fig. EB Even, when in this 
case the scaling for massless electrons is better than the corresponding one for Schrodinger 
electrons, see Figs. |S(d) and^d), the difference between the reference and the scaled curves 
is significant, compare the solid-black and the dotted-blue curves in Fig. E2(d). 

Even more interesting, the three scaling rules can be unified by means of the following 
general scaling expression: 

T(E, N, Do, w) « [T(—E, N + m, V 0 , -w)} 3 ^ 2 , (14) 

a k a 

where the factors k and a come in multiples of two, and m is the difference between gen¬ 
erations. As we can see in Fig. [TQI this general rule works pretty well for massive Dirac 
electrons. 

After analysing the different scaling rules, we obtain that, in general, the transmission 
properties of massive Dirac electrons present scalability, and as a consequence self-similarity. 
However, so far our comparison between the reference and scaled curves has been visual. 
What we need now is to get a quantitative measure of how similar are these curves. In order 
to address this issue, we compute the root mean square deviation or simply the rmsd, which 
is a standard measure to compare curves. If the curves tend to be exactly the same, this value 
is practically zero. For instance, the rmsd of massive Dirac electrons at normal incidence for 
the scaling between generations N, heights Vo and lengths w are 0.134303480, 0.110821329 
and 0.133542433, respectively. In the same way, the corresponding ones for Schrodinger 
electrons are: 0.256470233, 0.276356488 and 0.52796197. In the case of massive Dirac 
electrons at oblique incidence these values are 0.159395278, 0.131747648 and 0.157041878, 
while the corresponding ones for massless Dirac electrons become 0.448230058, 0.478004634 
and 0.463266015. Based on these results turns out that practically in all cases the rmsd is 
smaller for massive Dirac electrons. In particular, the rmsd for massive electrons is 3.6 times 
smaller than the rmsd of massless electrons in the case of the scaling between heights. This 
value becomes 2.5, if we compare massive Dirac electrons and Schrodinger electrons. At this 
point, it is also important to emphasize that if we change the energy range up to 1 eV, the 
matching between the reference and scaled transmission for massive Dirac electrons is even 
better. Specifically, the rmsd becomes one order of magnitude smaller. As a final remark, 
we can say that from both the qualitative and quantitative standpoints, the transmission 
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properties of massive Dirac electrons are self-similar. 


IV. CONCLUSIONS 

In summary, we have studied the self-similarity and scalability of the transmission proper¬ 
ties of Cantor Graphene Structures. The transfer matrix formalism has been used to compute 
the transmission spectra. In particular, three mathematical expressions are proposed and 
implemented to this end, finding three different scaling rules, which are associate with three 
fundamental parameters: the generation number, the height of the barriers and the length 
of the system. Thus, a comparative analysis of the scaling between generations, heights 
and lengths of the system for massive and massless Dirac electrons as well as Schrodinger 
electrons is carried out, obtaining that the better scaling is for massive Dirac electrons. 
This result is valid for both normal and oblique incidence. Finally, as far as we can see 
to obtain self-similar characteristics in the physical properties, such as the transmittance, 
are necessary three conditions: relativistic Dirac electrons, a self-similar structure and the 
non-conservation of the pseudo-spin (massive Dirac electrons). 
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FIG. 1. (Color online) Graphical representation of Substrate-Based Cantor Graphene Structures 
(SCGSs). (a) Representation of the cross-section of SCGSs. This structure is generated by the 
interaction between the substrate and the graphene layer. The graphene sheet is deposited above 
alternating substrates, such as Si02 (gray slabs) and SiC (blue slabs), obtaining regions with 
linear and parabolic dispersion relations, respectively, (b) Distribution of the Dirac cones and 
Dirac paraboloids along the structure. The bandgap in the Dirac paraboloids is E g = 2 1', the 
blue filled area indicates occupied states and the blue sphere describes an impinging massive Dirac 
electron at the Fermi energy, (c) Band-edge profile of SCGSs, the barriers have a height Vq = t 1 , 
and the widths of barriers and wells are dB 2 = w/9, dW 2 = w/9 and dWi = w/ 3. The cyan and 
blue lines correspond to electrons and holes, respectively. All the illustrations are presented for the 
generation N = 3. 
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FIG. 2. (Color online) Graphical representation of Electrostatic Cantor Graphene Structures 
(ECGSs). (a) Representation of the cross-section of ECGSs. This structure is generated when 
a graphene sheet is sited on a non-interacting substrate such as Si 02 (gray slab) along with a back 
gate (black slab) and alternating top gates (yellow slabs), (b) Dirac-cone distribution along the 
structure. In this case, Vo is proportional to the strength of the voltage applied to the back and 
top gates. The blue sphere describes an impinging massless Dirac electron at the Fermi energy, 
(c) Band-edge profile of ECGSs with Vo as height of the barriers. The cyan line corresponds to 
electrons, since we are considering positive voltages. All the illustrations are presented for the 
generation number N = 3. 
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FIG. 3. (Color online) Graphical representation of typical Al^Gai-^As/GaAs Semiconductor 
Structures (SSs). (a) Representation of the cross-section of SSs. This structure can be created 
by alternating semiconductors with dissimilar bandgaps, such as GaAs (gray slabs) and AlGaAs 
(blue slabs), (b) Dirac-paraboloid distribution, the blue filled area indicates occupied states and 
the blue sphere describes an impinging Schrodinger electron at the Fermi energy, (c) Band-edge 
profile of SSs. The cyan and blue lines correspond to electrons and holes, respectively. The height 
of the barriers is modulated by the aluminum concentration (x). All the illustrations are presented 
for the generation number N = 3. 
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FIG. 4. (Color online) Scaling between generations at normal incidence for SCGSs (first row) and 
SSs (second row), (a) Comparison of the transmittance as function of energy for generations N = 6 
(solid-black lines) and N = 7 (dashed-red lines), (b) Using the appropriate scaling equation (eq. 
ED, we can corroborate that the reference (N = 6) and scaled (N = 7) curves are approximately 
the same, solid-black and dotted-blue lines, respectively, (c) and (d) correspond to the same 
comparison as in (a) and (b) but for SSs, or equivalently Schrodinger electrons. In this case the 
scalability is not present. The structural parameters of both systems SCGSs and SSs are w = 1000a 
and Vq = 0.13 eV. 
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FIG. 5. (Color online) Scaling between generations at oblique incidence (6 = 7r/4) for SCGSs (first 
row) and ECGSs (second row), (a) Comparison of the transmittance as function of energy for 
generations N = 6 (solid-black lines) and N = 7 (dashed-red lines), (b) Using the appropriate 
scaling equation (eq. flTT) . we can corroborate that the reference (N = 6) and scaled (N = 7) curves 
are approximately the same, solid-black and dotted-blue lines, respectively, (c) and (d) correspond 
to the same comparison as in (a) and (b) but for ECGSs, or equivalently massless Dirac electrons. 
In this case the scaling properties are not visible. The structural parameters are the same as in 
Fig. m 
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FIG. 6. (Color online) Scaling between heights of the barriers at normal incidence for SCGSs (first 
row) and SSs (second row), (a) Comparison of the transmittance as function of energy for the 
heights Vo = 0.10 eV (solid-black lines) and Vo = 0.05 eV (dashed-red lines), (b) By implementing 
the appropriate scaling equation (eq. Il2l) . we can see that the scaled curve corresponding to 
Vo = 0.05 eV is practically similar to the reference one Vo = 0.10 eV, dotted-blue and solid-black 
lines, respectively, (c) and (d) correspond to the same comparison as in (a) and (b) but for SSs. 
In this case the scalability is not present. The structural parameters for both systems, SCGSs and 
SSs, are w = 1000a and N = 6. 
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FIG. 7. (Color online) The same as in Fig. [H but for SCGSs (first row) and ECGSs (second 
row), at oblique incidence 8 = it /A. As in Fig. [6] massive Dirac electrons present excellent scaling 
properties. 
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FIG. 8. (Color online) Scaling between lengths of the system at normal incidence for SCGSs (first 
row) and SSs (second row), (a) Comparison of the transmittance as function of energy for lengths 
w = 1000a (solid-black lines) and w = 500a (dashed-red lines), (b) By applying eq. [T3] the scaled 
curve (w = 500a) is obtained, which has an excellent match with the reference one (w = 1000a). 
(c) and (d) correspond to the same comparison as in (a) and (b) but for SSs. In this case the 
scalability is not present. Here, the structural parameters are N = 6 and Vq = 0.13 eV. 
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FIG. 9. (Color online) The same as in Fig. 0 but for SCGSs (first row) and ECGSs (second 
row), at oblique incidence 8 = 7t/4. As in Fig. 0 massive Dirac electrons present excellent scaling 
properties. 
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FIG. 10. (Color online) General scaling at normal incidence for SCGSs. (a) Comparison of the 
transmittance as function of energy for N = 6, Vo = 0.10 eV and w = 1000a (solid-black lines) and 
N = 7, Vo = 0.05 eV and w = 500a (dashed-red lines), (b) By applying eq. HHthe scaled curve 
(dotted-blue lines) is obtained, which has an excellent match with the reference one (solid-black 
lines). As in all previous figures massive Dirac electrons present outstanding scaling properties. 
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